A New Family of Light Beams and Mirror Shapes for Future LIGO Interferometers 
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Advanced LIGO's present baseline design uses arm cavities with Gaussian light beams supported 
by spherical mirrors. Because Gaussian beams have large intensity gradients in regions of high 
intensity, they average poorly over fluctuating bumps and valleys on the mirror surfaces, caused 
by random thermal fluctuations (thermoelastic noise). Flat-topped light beams (mesa beams) are 
being considered as an alternative because they average over the thermoelastic fluctuations much 
more effectively. However, the proposed mesa beams are supported by nearly fiat mirrors, which 
experience a very serious tilt instability. In this paper we propose an alternative configuration in 
which mesa-shaped beams are supported by nearly concentric spheres, which experience only a 
weak tilt instability. The tilt instability is analyzed for these mirrors in a companion paper by 
Savov and Vyatchanin. We also propose a one-parameter family of light beams and mirrors in 
which, as the parameter a varies continuously from to ir, the beams and supporting mirrors get 
deformed continuously from the nearly fiat- mirrored mesa configuration ("FM") at a — 0, to the 
nearly concentric-mirrored mesa configuration ("CM") at a = it. The FM and CM configurations at 
the endpoints are close to optically unstable, and as a moves away from or n, the optical stability 
improves. 

PACS numbers: Valid PACS appear here 
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INTRODUCTION 



The initial gravitational-wave detectors in the Laser 
Interferometric Gravitational wave Observatory (LIGO) 
are now near design sensitivity and are taking science 
data [1]. The interferometers will be upgraded to a much 
more sensitive advanced-LIGO design beginning in about 
2007. Until 2003, the baseline design for advanced LIGO 
used nearly flat but spherical mirrors in its arm cavities. 
However, in 2003, Sidles and Sigg [2] showed that these 
mirrors experience a strong tilt instability: when the mir- 
rors are tilted symmetrically, the light beam slides across 
their surfaces to an off-center location and its light pres- 
sure then pushes hard to increase the tilt. Sidles and Sigg 
proposed switching to mirrors that arc segments of nearly 
concentric spheres (radii of curvature slightly larger than 
half the cavity length); such mirrors, they showed, can 
support Gaussian beams of the same (large) radius as 
the baseline design, while experiencing a much weakened 
tilt instability. This triggered a change of the baseline 
design for advanced LIGO to nearly concentric, spherical 
mirrors. 

Gaussian beams have the serious disadvantage that, 
because of their steep intensity gradient over most of 
the beam's area, they average poorly over the fluctu- 
ating bumps and valleys on the mirrors' surfaces that 
are caused by thermal fluctuations (thermoelastic noise). 
O'Shaughnessy and Thorne [3-5] have proposed improv- 
ing the averaging and thereby reducing the thermoe- 
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lastic noise substantially, by replacing the arm cavities' 
Gaussian beams by flat-topped beams (mesa beams, as 
Willems has named them), which are supported by nearly 
flat, Mexican-hat-shaped mirrors. O'Shaughnessy, Stri- 
gin and Vyatchanin [5, 6] have shown that the thermoe- 
lastic noise for mesa beams is three times weaker in noise 
power than for the baseline Gaussian beams (with beam 
sizes that produce the same diffraction losses), and cor- 
respondingly that the mesa-beam interferometers could 
see farther into the universe, producing event rates for 
inspiraling binaries about three times higher than the 
baseline Gaussian-beam design. D'Ambrosio et. al. [5-7] 
have shown that the nearly flat mesa-beamed mirrors are 
practical in all respects that could be analyzed theoreti- 
cally, and D'Ambrosio, DeSalvo, Simoni and Willems [8] 
are building a prototype optical cavity with mesa beams, 
to explore practical issues experimentally. 

Unfortunately, the nearly flat, Mexican-hat-shaped 
mirrors ("FM") proposed by O'Shaugnessy and Thorne 
to support mesa beams (Sec. II below), like the nearly 
flat, spherical mirrors of the pre-2003 baseline design, ex- 
perience a severe tilt instability (Vyatchanin [9]; Savov 
and Vyatchanin [10]). In this paper, motivated by the 
Sidles-Sigg result that, for Gaussian beams and spher- 
ical mirrors, the tilt instability is greatly weakened by 
switching from nearly flat to nearly concentric mirrors, 
we propose (Sec. Ill) a new, nearly concentric mirror 
design ("CM") that supports mesa beams. In a com- 
panion paper, Savov and Vyatchanin [10] show that the 
tilt instability is weaker for these CM mirrors than for 
any other mirrors thus far considered — FM, nearly-flat 
spherical, and nearly-concentric spherical. 

In Sees. II and III, we mathematically construct our 
FM and CM beams and their Mexican-hat mirror shapes 
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FIG. 1: Optical axes of the families of minimal Gaussians 
beams used to construct: (a) an FM mesa beam [6], denoted 
in this paper a = 0; (c) our new CM mesa beam, denoted 
a = 7r; (b) our new family of hyperboloidal beams, which 
deform, as a varies from to it, from a FM beam (a) into a 
CM beam (c). 



by superposing minimal-radius Gaussian beams with op- 
tic axes that are the generators of cylinders (for FM) 
and of cones (for CM); Fig. la,c. In Sec. IV we introduce 
a one-parameter family of "hyperboloidal" light beams 
and supporting mirrors, computed by superposing min- 
imal Gaussians whose optic axes are the generators of 
hyperboloids; Fig. lb. For each hyperboloidal beam, the 
hyperboloid's generators (minimal-Gaussian optic axes) 
have a fixed twist angle a. As a is varied continually, the 
hyperboloidal light beams deform continually from mesa- 
shaped FM form (at a = 0, where the hyperboloids are 
cylinders) to sharply peaked Gaussian form (at a = n/2), 
to mesa-shaped CM form (at a = n, where the hyper- 
boloids degenerate to cones). 



[4, 6]: One superposes minimal Gaussian beams 1 with 
their optic axes all parallel to the cavity axis and dis- 
tributed uniformly over a disk with some radius D, as 
shown in Fig. la. Each minimal Gaussian beam (field) is 
given by 
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Here w is transverse distance from the beam's optic axis; 
Q is distance parallel to the optic axis with Q = at the 
beam waist; £ = L/2 is half the length of LIGO's arm 
cavity (2 km) and is also equal to the beam's Rayleigh 
range; b = y/\L/2ir = y/\£/ir = 2.603 cm (with A = 
1.064/im the light wavelength) is the radius, at the 1/e 
point of the beam's intensity distribution, at the ends of 
the cavity, i.e. at £ = £; and b is also the radius, at the 1/e 
point of the beam's amplitude distribution, at the beam's 
waist, C = 0. Note that the last phase factor, 2££/b 2 , is 
actually k( in disguise, with k = 2n/X the light's wave 
number. We adjust A or £ slightly so that at £ = £ and 
zu = 0, is real and positive, i.e. the sum of the last two 
phase factors is a multiple of 27r. Then in the immediate 
vicinity of the mirror plane, at £ = £ + 5Q (with \8C\ <C 6), 
the minimal Gaussian has the simple form 
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with k = 2£/b 2 . (Here and throughout this paper we ig- 
nore fractional corrections of order A/6 <~ b/£ ~ 10 -5 .) 
The mesa-beam (FM-beam) field, constructed by super- 
posing minimal Gaussians as in Fig. la, is given by (Sec. 
IIA of [5]) 



U (r,z,D) 



*(%/ (x - x ) 2 + (y - y ) 2 , z) dx dy 



(3) 

Here r = a/x 2 + y 2 is radius from the cavity's central 
axis, the integral is over Cartesian coordinates (x , y ) of 
the Gaussians' optic axes, and the integral extends over 
the interior of the diskC^i with radius D, i.e. \J x 2 + y 2 < 
D. The subscript on Ua is the value a = of the twist 
angle of the Gaussians' optic axes, when one regards this 
FM beam from the viewpoint of the hyperboloidal family 
of beams (Sec. IV below). 

By inserting expression (2) with 6( = into Eq. (3), 



II. MESA BEAMS SUPPORTED BY NEARLY 
FLAT MIRRORS (FM BEAMS; a = 0) 

The mesa beams supported by nearly flat Mexican-Hat 
mirrors ("FM" beams) can be constructed mathemati- 
cally by a procedure due to O'Shaughnessy and Thorne 



1 By "minimal Gaussian beam" we mean the fundamental, TEM00 
mode of an optical resonator with spherical mirrors, with the 
mirror radii of curvature adjusted so the Gaussian-shaped in- 
tensity distributions on the two mirrors are identical and have 
the minimum possible radii b at the 1/e point of the intensity 
distribution. 



3 



/a = 0.87t \ 




FIG. 2: The correction H a (r) to the mirror shape for hy- 
perboloidal beams with D = 10 cm and with twist angles a 
between 7r/2 and n. For a — 0, the correction is the negative 
of that for a = n; for a — O.ln it is the negative of that 
for 0.97r; for any a between and tt/2, it is the negative of 
that for tv — a. For a — tt (the Mexican-hat correction for 
our new CM mesa beam), Ho(r) drops to about —500 nm 
(half the wavelength of the light beam) at r = 16 cm (the 
mirror's edge). These corrections are added onto the fiducial 
spheroidal shape S a (r) [Eq. (13)]. 



we obtain for the FM beam at the mirror plane z = I , 

U (r,£,D) = (4) 



exp 



-[{x-x ) 2 + {y-y ) 2 ][l-i] 



2b 2 



dx dy 



The mirror surface must coincide with a phase front 
of this mesa beam (FM beam), i.e. it must have a shape 
z = £+H (r) such that axg[Uo(r, £+H , D)] = constant = 
axg[Uo(0,£, D)]. In the vicinity of the mirror the phase 
of each minimal Gaussian varies as k5Q = k5z [Eq. (2)], 
so axg[U (r,£ + H ,D)] = arg[t/ (r, £, D)] + kH , and the 
shape of the mirror surface must be 



Ho(r) 



{axg[U a (0,i,D)-axg\U (r,l,D)]} . (5) 



By carrying out the integral (4) analytically in one di- 
mension then numerically in the other, and then insert- 
ing into Eq. (5), O'Shaughnessy and Thorne find the 
"Mexican-hat" mirror shape H(r) shown as a solid line 
in Fig. 2. 

To high accuracy, the field Uq on the mirror surface 
differs from that on the plane z — £ only by the phase 
factor e tkH o( r ) ; so th e intensity distribution on the mirror 
is the same as at z — b, i.e., it is 



J (r) cx \U (r,£,D)f 



(6) 



This intensity has the mesa shape shown as a solid line 
in Fig. 3. 



III. MESA BEAMS SUPPORTED BY NEARLY 
CONCENTRIC MIRRORS (CM BEAMS; a = tt) 

Our proposed mesa beams with nearly concentric, 
spherical mirrors ("CM" beams) can be constructed by 
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FIG. 3: The light beam's un- normalized intensity \U a \ as a 
function of radius r on the mirror, for hyperboloidal beams 
with D — 10 cm and various twist angles a. For a — and 7r, 
the intensity has the mesa shape; for a = 0.5 it is a minimal 
Gaussian. 



overlapping minimal Gaussians whose optic axes all pass 
through the center of the cavity [Fig. 1(c)], and are 
distributed uniformly inside a cone with angular radius 
G = D/£ 

It should be clear from this construction that the re- 
sulting mesa beam will have a beam radius approximately 
equal to D at the mirrors (z = £), and approximately 
equal to b (the minimal Gaussian radius) at the cavity's 
center, z = 0. For advanced LIGO, D is approximately 
4& [5], so the waist of this mesa beam will be approxi- 
mately 4 times narrower than the beam on the mirrors. 
This contrasts with a mesa beam supported by nearly 
flat mirrors (previous section) , for which the waist is only 
slightly narrower than the beam on the mirrors. 

Near the mirrors, the phase fronts of this CM beam 
will be nearly concentric spheres centered on the point 
(r, z) — (0, 0) through which the Gaussians' optic axes 
pass, so we shall evaluate the CM field as a function of 
radius r on this fiducial sphere, 



z = S„(r) = \/£ 2 - r 2 ~ £ - r 2 /2£ 



(7) 



[Here and below we use a subscript n to denote mesa- 
beam quantities with nearly concentric mirrors, i.e. CM 
quantities. This is because the minimal Gaussians used 
to generate the CM beam have twist angles a = tt; 
see Sec. IV.] For each minimal Gaussian, this fiducial 
sphere bends away from the Gaussian's transverse plane 
by an amount 6£ — — / cu 2 /2£, so on this fiducial sphere 
the Gaussian's phase factor k5( — (2£/b 2 )S( is equal to 
—vj 2 /b 2 . As a result, the Gaussian field on the fiducial 
sphere is [cf. Eq. (2)] 
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Correspondingly, these minimal Gaussians superpose, on 
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the fiducial sphere, to produce a CM field given by 
U„(r,S v ,D)= (9) 



L 



exp 



.[( X - Xo f + (y-y f][l + t ] 

2b 2 



dx dy . 



Notice that this [^(r, S„, D) is the complex conjugate 
of the FM field Uo(r,£,D) evaluated on the transverse 
plane z — I (the fiducial surface for the case of nearly 
flat mirrors); Eq. (4). 

As for the nearly flat (FM) case, the phase of the 
CM field will vary, with distance 5z from the fiducial 
sphere, nearly proportionally to k5z; and correspond- 
ingly, the mirror's surface, Sz — H n (r) (a surface of con- 
stant phase), will be given by the analog of Eq. (5): 

HAr) = k- 1 {arg[[^(0, S n , D) - arg[*7 w (r, S„,D)}} . 

(10) 

Because U^{vj, S n , D) is the complex conjugate of 
U (w,£,D), Eqs. (5) and (iO) imply that 



H„(r) = -H (r) . 



(11) 



In words: to support mesa beams with the same beam 
size D on their mirrors, the nearly concentric mirrors 
and the nearly flat mirrors must deviate from precisely 
concentric spheres z = S^ir) and precisely flat planes 
z = £ by equal and opposite displacements Sz = H(r). 
This fact was discovered in numerical work by one of 
us (MB) and was later proved numerically in a much 
wider context by Savov [iO] and analytically by Agresti, 
d'Ambrosio, Chen and Savov [11], before we found the 
above demonstration. 

Because [to the accuracy of our analysis, 0(X/b)] the 
field U n is the same, aside from phase, on the mirror 
surface as on the fiducial sphere S n , the light's intensity 
distribution is the same on the mirror as on S-n- 



(12) 



Moreover, because C/ W (r, Sn,D) is the complex conju- 
gate of Uo(r,£,D), they have the same moduli and in- 
tensity distributions — i.e., the CM beam has the same 
mesa-shaped intensity distribution as the FM beam (solid 
curve in Fig. 3 below). This fact was discovered in nu- 
merical work by one of us (MB) and was later proved nu- 
merically in a much wider context by Savov [10] and an- 
alytically by Agresti, d'Ambrosio, Chen and Savov [11], 
before we found the above demonstration. 



IV. HYPERBOLOIDAL BEAMS SUPPORTED 
BY NEARLY SPHEROIDAL MIRRORS 

One can smoothly transform the FM beams into CM 
beams, and the in-between beams may be interesting for 
L1GO. In this section, we will focus on one way to make 
such a transformation. 

We will first look at a smooth deformation of the ge- 
ometric body formed by the optic axes of the minimal 



Gaussians that are used in constructing the FM and CM 
beams. For a FM beam, the axes of the minimal Gaus- 
sians lie on coaxial cylinders, while for CM beams they lie 
on coaxial cones. It is well-known that one can smoothly 
deform a cylinder into a cone as follows. The generators 
of a cylinder of height 2£ and radius r (Fig. la) arc lines 
that join points with cylindrical coordinates (r, <f), —£) on 
the base circle to points (r, <f), £) on the top circle. The 
generators of a symmetric cone of height 2£ and end radii 
r (Fig. lc) are lines that join points (r, <j>, —£) and points 
(r,(f>+n,£). A path from the cylinder to the cone is given 
by a family of hypcrboloids generated by lines that join 
points (r,<f>,—£) and points (r, <j> + ot,t) (Fig. lb). For 
a = one obviously gets the cylinder and for a = n, the 
cone. 

We therefore propose constructing a new two- 
parameter family of light beams, and the mirrors that 
support these beams, using the O'Shaugnessy-Thorne 
technique of superposing minimal Gaussians. The pa- 
rameters are {a, D}, and for given values of {a,D} the 
minimal Gaussians have their optic axes uniformly dis- 
tributed on the hyperboloid generators that reach from 
(r, </>, —£) to (r, (f> + a, £) (Fig. lc), with <f> running from 
to 2tt and r confined to the interior of the disk Co, r < D. 
For a = these hyperboloidal beams will be mesa beams 
with nearly flat mirrors, i.e. FM beams. For a = n, they 
will be mesa beams with nearly concentric mirrors, i.e. 
CM beams. 

We can construct explicit expressions for the shapes 
of the mirrors that support these hyperboloidal beams, 
and expressions for the fields on those mirrors, using the 
same method as in the FM case (Sec. II) and the CM case 
(Sec. Ill): Because the phase of each minimal Gaussian 
varies nearly proportionally to the surface of constant 
phase at the mirror location will be nearly the same as the 
"fiducial" surface obtained by cutting off each Gaussian's 
optic axis at ( — £. One can show that, with the optic 
axes being generators of hyperboloids, the surface formed 
by their ends at constant distance £ = I from the cavity's 
mid point is the fiducial spheroid 

z = S a (r) ee y/p r 2 sin 2 ( Q /2) ~ £ - ^ ,2) . 

(13) 

We can compute our hyperboloidal field U a (r, S a , D) 
on this fiducial spheroid by superposing our minimal 
Gaussians with the aid of Fig. 4. In this figure P is the 
point on the spheroid S a at which we wish to compute the 
field. The vector p reaching from the spheroid's center 
point O (the center of our hyperboloidal field's cross sec- 
tion) to P has Cartesian coordinates p = (r, 0, Z), where 
Z = S a (r) - £ = -(r 2 /2£)sm 2 (a/2). The optic axis of 
a minimal Gaussian, over which we will integrate, inter- 
sects S a at the point Q, which has Cartesian coordinates 
q = (r o cos0 o ,r o sin0 o ,Z o ), where Z a = S a (r a ) - £ = 
-(r 2 j2£)sm 2 (a/2). The optic axis of this minimal Gaus- 
sian points along the unit vector n = {(r /2f)[cos(</> ) — 
cos(0 o — a)], {r /2£)[sa\(<p ) — sin(</> — a)], 1}. [Here as 
elsewhere we neglect corrections of order r/£ ~ r a /£ <~ 
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FIG. 4: Geometric construction for computing the hyper- 
boloidal field U a (r, S a , D) on the fiducial spheroid S a (a seg- 
ment of which is shown dotted). 



bjl ~ A/6 ~ 10~ 5 .] The vector s = p — (q+n<5£) reaches 
orthogonally from the minimal Gaussian's optic axis to 
the point P. The length of this vector is the radius w of 
P as measured in cylindrical coordinates centered on the 
minimal Gaussian's optic axis, 



|p - q| = \A 



2rr D cos < 



(14) 



where the second expression, accurate to O(Xfb), can be 
deduced from the above equations. The distance 6( along 
the optic axis n, at which the normal s intersects the axis, 
is determined by the orthogonality relation s • n = 0: 



-1 



w 2 sin 



rr n sin a sin < 



(15) 

The field U a (r, S a , D) on the spheroid S a is obtained 
by adding up the minimal Gaussians (2) with w and 5£ 
given by Eqs. (14) and (15), and with k = 2£/b 2 , and by 
then doing some simple algebra: 







dr D d 


</> exp 


Jo Jo 





U a (r,S a ,D) = 

(r 2 + r 2 - 2rr a cos0 o ) 
Zb 2 



b 2 

(1 — i cos a) 



i-TTT sin <po sma 



(16) 



The radial integral can be carried out analytically yield- 
ing an expression involving error functions, and the an- 
gular integral can then be done numerically. 

The field (16) cannot be sensitive to the chirality of 
the optic axes' twist, i.e. to the sign of a, since it is a 
scalar complex function of r: U- a = U a . This tells us 
that the relevant range for a is to tt. Replacing a by 
7T — a and changing the sign of a is equivalent to complex 
conjugating U a ; therefore: 



(17) 



For a — 0, the fiducial spheroid So(r) is the transverse 
plane and the field (16) is the FM mesa beam U$ [Eq. 
(4)]. For a = ir/2, the fiducial spheroid SV/2 is a sphere 



of radius R = L = 2£ (the distance between the mir- 
rors), and both the radial and the angular integrals can 
be carried out analytically, giving for the field on that 
sphere 



U- 



zr/2 



constant exp [— r 2 /2& 2 ] 



(18) 



this is precisely the minimal Gaussian beam [Eq. (2) with 
w = r, evaluated at kS( = -kr 2 /2R = -(2£/b 2 )r 2 /U = 
—r 2 /2b 2 ]. For a — ir, the fiducial spheroid S n (r) is a 
sphere with radius I = L/2 and the field (16) is the CM 
mesa beam [Eq. (9)]. Thus, as a varies from to 7r, 
U a deforms continuously from the FM mesa beam a = 
0, through a set of hyperboloidal beams to a minimal 
Gaussian at a — 7r/2, and on through another set of 
hyperboloidal beams to the CM mesa beam a = ir. 

As for the FM and CM beams, so also for the hyper- 
boloidal beam (16) (and for the same reasons), the mir- 
ror's surface must be displaced longitudinally from the 
fiducial spheroid z = S a (r) by Sz — H a (r), where 

H a (r) = k- 1 {arg[C/ Q (0, S a ,D)- arg[C/ Q (r, S a ,£>)]}. 

(19) 

This equation and Un-a = —U^ Pq- (17)] tell us that 

H v - a {r) = -H a (r) . (20) 

This is a special case of a duality relation discovered nu- 
merically by Savov and Vyatchanin [10] and proved an- 
alytically by Agresti, d'Ambrosio, Chen and Savov [11]. 
Figure 2 shows these mirror shape corrections for various 
a's. The light intensity on the mirrors is given by the 
obvious analog of Eqs. (6) and (12): 

I a (r) oc \U a (r,S a ,D)\ 2 . (21) 

This equation and C/^_ a = — U* [Eq. (17)] tell us that 

/*-»(»•)= ia(r) . (22) 

This is another special case of the duality relation discov- 
ered by Savov and Vyatchanin [10], and proved analyti- 
cally by Agresti, d'Ambrosio, Chen and Savov [11]. The 
intensity distribution (22), (16) is shown, for various a, 
in Fig. 3. 

V. CONCLUSIONS 

For twist angles a near and 7r, the hyperboloidal 
beams introduced in this paper have the flat-top form 
needed to reduce thermoelastic noise in LIGO. The ra- 
dius of the flat top is largest for a — and a — ir (the 
FM and CM mesa beams) and smallest for a — it/ 2 (the 
minimal Gaussian). 

Because the mirrors are most nearly flat or concentric 
for the mesa configurations, a = or 7r, those configura- 
tions are most nearly optically unstable. (Near instabil- 
ity goes hand in hand with large beams on the mirrors, 
which are needed to control thermoelastic noise.) 
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The results of Savov and Vyatchanin [11] suggest that 
the tilt instability is smallest for a = 1 and worst for 
a = ir. 

These considerations suggest that the optimal config- 
uration for advanced LIGO will be near a = n, but 
whether the optimum is precisely at a — n (the CM con- 
figuration) or at some modestly smaller a will depend 
on practical and thermoelastic-noise considerations not 
examined in this paper. 
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